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ABSTRACT 

We present the latest results from a spectroscopic survey designed to uncover the 
i-C ■ hidden population of AM Canum Venaticorum (AM CVn) binaries in the photometric 

database of the Sloan Digital Sky Survey (SDSS). We selected ^2000 candidates based 
on their photometric colours, a relatively small sample which is expected to contain 
the majority of all AM CVn binaries in the SDSS (expected to be ^50). 

We present two new candidate AM CVn binaries discovered using this strategy: 
SDSS J104325.08+563258.1 and SDSS J173047.59+554518.5. We also present spec- 
tra of 29 new cataclysmic variables, 23 DQ white dwarfs and 21 DZ white dwarfs 
discovered in this survey. 

The survey is now approximately 70 per cent complete, and the discovery of 
seven new AM CVn binaries indicates a lower space density than previously predicted. 
From the essentially complete g ^ 19 sample, we derive an observed space density of 
(5 ± 3) x 10~ 7 pc~ 3 ; this is lower than previous estimates by a factor of 3. 

The sample has been cross-matched with the GALEX All-Sky Imaging Survey 
database, and with Data Release 9 of the UKIRT (United Kingdom Infrared Telescope) 
Infrared Deep Sky Survey (UKIDSS). The addition of UV photometry allows new 
colour cuts to be applied, reducing the size of our sample to ~1100 objects. Optimising 
our followup should allow us to uncover the remaining AM CVn binaries present in the 
SDSS, providing the larger homogeneous sample required to more reliably estimate 
their space density. 

Key words: accretion, accretion discs - binaries: close - stars: individual: SDSS 
J104325.08+563258.1, SDSS J173047.59+554518.5 - novae, cataclysmic variables - 
white dwarfs. 



1 INTRODUCTION 

The AM Canum Venaticorum (AM CVn) binaries are a 
class of ultracompact systems that consist of a white dwarf 
accreting helium rich material from a (semi-)degenerate 
donor. They are characterised by their short orbital peri- 
ods, which range from 5 to ~65 minutes, and an absence 
of hydrogen in their spectra. It is the degenerate nature 
of the mass donor that allows their periods to lie well be- 
low the believed minimum period spike (~80 min) of hydro- 
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gen rich cataclysmic variables (CVs: iRappaport et al.lll982l ; 
iGansicke et al.ll2009l ). Only 33 members of the rare AM CVn 
binary cl ass have been repor ted in the literature (the most 
recent bv lLevitan et al.ll2012l ) , 7 of which were discovered via 
systematic searches in the Sloan Digital Sky Survey (SPSS ; 
York et alJliuOol : iRoelofs et all 120051 : lAnderson etaD 120051 . 
20081). A recen t review of the AM CVn binaries is given by 
Solheiml (|201Ch . 

The observational features of an AM CVn binary are 
thought to depend on its orbital period, and it will pass 
through several distinct phases as the system evolves from 
period minimum towards longer orbital periods. In the 
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shortest period systems (P or b <~10 min) the accretion 
stream impacts dire ctly onto the surface o f the accretor, 
and no disc forms (|Marsh fc Steeghsl |2002| ; iRoelofs et al.l 
2010). For orbital periods of ~10 - ~20 minutes, the ac- 
cretion disc is in a stable 'high' state, and spectra are 
domi nated by helium absorp tion from the optically thick 
disc l|0'Donoghue et aH fl994). In the long period systems 
(forb >~40 min) the disc is in a stable low state; these sys- 
tems typically lack photometric variability, and th eir opti- 
cal spectra are dominate d by helium emission lines (|Warnerl 
Il995bl ; iRuiz et "al"1l200lh . The intermediate period systems 
(20 ^ Porb ^ 40 min) have unstable discs, and their ap- 
pearance varies between that of the high-state and the low- 
state systems, analogously to the hydrogen-r i ch dwarf novae 
(iTsugawa fc Osakil 1 19971 : iKotko et al] |2012| ; iRamsav et ail 
I2012T ). 

The AM CVn binaries are of great interest for 
many aspects of astrophysics. These systems represent 
the end product of seve r al finely-tuned evolut i onary 
pathways dlben fc Tutukovl Il99ll ; iNelemans et all l200ll ; 
iPodsiadlowski et aLlbOoU T and as such are of great inter- 
est for binary stellar evolution theory. They are also of 
particular interest due to their gravitational wave emis- 
sion, being some of the strongest known sources that would 
be detected by future low-frequency, space-borne gravita- 
tional wave detectors, the brighte st systems acting as cal- 
ibrators for such an experiment ( Stroccr & Vccchio |2006| ; 
IRoelofs et ai1l2006bl ; lNelemansll2009l ; iNissanke et al.l [201^ V 
Determining an accurate space density will allow improved 
constraints to be placed on the expected gravitational wave 
signal. The AM CVn binaries are closely related to the 
double degenerate pathway to Type la supernovae, and 
may themselves contribute to the supernova population 
(|Solheim fc Yungelsonll2005l ; iRuiter et al.ll2009h . They may 
also harbour dynamical time scale helium fus ion, producing 
rare s ub-luminous SN la -like events f'SN.Ia'; iBildsten et al.l 
2007; Bro wn et al.ll201lh . Establishing the space density is 
important for determining t he stability of mass transfer in 
double white dwarf binaries (|Marsh et al.ll2004l ). which are 
potential pro genitors for Typ e la supernovae, sdB stars and 
R CrB stars (|Webbinkl \l984 ). 

The first AM CVn binaries were discovered serendipi- 
tously in a number of different ways. INelemans et al.ll200ll 
conducted a population synthesis based on the two main 
proposed formation channels for AM CVn binaries, conclud- 
ing that improved observations were required in order to 
learn more about the AM CVn population. The discovery 
of six AM CVn binaries in the SDSS spectrosc opic database, 
via t h eir helium emission d o minat ed spectra (jRoelofs et al.l 
120051 ; lAnderson et al.1 120051 . |200ST ). provided the first rela- 
tively well defined, homogeneous sample, that allowed esti- 
mation of the local space density by calibrating the theo- 
retical popjilB ; tion_mo^els. This gave a value of 1-3 x 10 -6 
pc -3 (|Roelofs et al.l l^007ch . This was an order of magni- 
tude lower than the expected space density at the time, de- 
rived from populatio n synthesis (2 x 10 -5 pc -3 , based on 
INelemans et al.ll200lf ). but the accuracy was limited by the 
small sample size. 

Since the known AM CVn binaries have been found 
to occupy a relatively sparsely-populated region of colour 
space, which has a low spectroscopic completeness in the 
SDSS database, we have initiated a dedicated spectroscopic 
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Figure 1. Completeness of SDSS DR9 spectroscopic follow-up 
as a function of colour, u — g and g — r, to a limiting magnitude 
g = 20.5 (dereddened). The long period SDSS AM CVn binaries 
are indicated by star symbols. The solid line marks the blackbody 
cooling track, the dotted and dot-dashed lines indicate model 
cooling sequences for DA and DB white dwarfs. 



survey designed to uncover the 'hidden' population of AM 
CVn binaries in the SDSS photometry. Details of this pro- 
gramme, including the sample selection crit eria, and the 
first n ew AM CVn binary, were presented by IRoelofs et al. | 
2009); four further discoveries were described bv lRau et al.l 



2010) 



The sample has since been extended to include new tar- 
gets from SDSS Data Release 7. This increase in the survey 
area should lead to a corresponding increase in the size of the 
resulting AM CVn sample (which is desirable given the still 
small number of systems found), however, it also increases 
the observation time required to complete the programme. 
We therefore investigate the possibility of increasing our AM 
CVn detection efficiency using new colour-colour cuts based 
on photometry at longer and shorter wavelengths than the 
five SDSS bands provide. Here we present an update on the 
status of the programme, detail our findings so far, and pro- 
pose improved selection criteria to increase our efficiency. 

We describe the status of the survey and our observa- 
tions in Section 2. In Section 3 we discuss cross-matching 
to the GALEX All-Sky Imaging Survey and the UKIRT 
(United Kingdom Infrared Telescope) Infrared Deep Sky 
Survey. In Section 4 we discuss our identification procedure, 
and present our CV, DQ white dwarf and DZ white dwarf 
samples. In Section 5 we detail our results, present two new 
candidate AM CVn binaries, describe new color cuts based 
on GALEX photometry, and report our revised value for the 
space density. We discuss these results in Section 6. 



2 SURVEY OBSERVATIONS 

Fig. [T] shows the spectroscopic completeness of the SDSS 
Data Release 9 as a function of u — g and g — r 
colour, to a limiting (/-band magnitude of 20.5. Also 
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plotted is the blackbody cooling track, and Bergeron 
model colours of the cooling sequences of hydrogen (DA) 

and helium (DB) atmospher e white dwarfs (log g = 

8.0; iHolberg fc Bergeron! |2006|; iKowalski fc Saumonl [20061 ; 
iTremblav et al.ll201ll ; lBergeron et al.ll201lh . The known AM 
CVn binaries lie in an area of relatively low spectroscopic 
coverage, largely offset from the DB white dwarf track. The 
areas to the lower left and right are targeted for DA white 
dwarfs and quasars^ 

As shown by iRoelofs et ail (|2009h . AM CVn binaries 
with a He I 5875 equivalent width ^30 A are easily de- 
tectable at a signal-to-noise ratio of 10, and resolutions as 
low as R = 300. Even systems with weaker lines should 
still show some sign of helium emission that would be con- 
firmed by a second, higher quality spectrum. The possibility 
that there e xist some syst e ms w ith much weaker lines was 
discussed bv lRoelofs et al.l (|2009T ). They concluded that the 
majority of the population should be detectable. 

Low-resolution, low signal-to-noise ratio spectra of ob- 
jects in our colour-selected sample have been obtained with 
several different instrument setups, using a variety of low 
and intermediate resolution spectrographs. To date, spec- 
troscopic observations of 1403 of the 1947 candidates have 
been completed (see Fig. These data were obtained us- 
ing a number of telescopes: the Isaac Newton Telescope 
(INT), William Herschel Telescope (WHT), and Nordic Op- 
tical Telescope (NOT), on the island of La Palma; the 1.5- 
m Tillinghast telescope at the Fred Lawrence Whipple Ob- 
servatory, Mt. Hopkins, Arizona; the Very Large Telescope 
(VLT), at Paranal; the 200- inch Hale telescope at the Palo- 
mar Observatory; both Gemini telescopes, South in Chile 
and North on Hawaii; and the New Technology Telescope 
(NTT), at La Silla. The largest telescopes, VLT and Gem- 
ini, have been used to observe the faintest part of the sample, 
(/-band magnitude > 20. A log of our observations, listing 
the instrument setups used and the numbers of candidates 
observed, is given in Table Q] 

Instrument setups were chosen such that all spectra 
cover the 4000-7000 A wavelength range that includes the 
predominant lines of both helium and hydrogen, allowing 
identification of AM CVn binaries from their strong helium 
emission, and lack of hydrogen. 

The spectra obtained with the FAST spectrograph were 
reduced using the spectral extraction pipeline provided by 
the observato ry. This pipelin e is ba sed on standard IRAF 
routines, see iTokarz fc Roil (| 19971 ). All other data were 
reduced using optimal extra ction as implemented in the 
Pamel^I code l|Marshl Il989), and the Starlink packages 
Kappa, Figaro and Convert. Wavelength calibration was 
obtained from various arc lamp exposures taken each night. 
Flux calibration was achieved with various standard stars 
observed at the beginning or end of each night. The flux 
calibration is not absolute, and approximately one third of 
the spectra do not have flux calibration. 

Fig. [2] shows the current spectroscopic completeness of 
our sample. The target list at the bright end of the distri- 
bution has been effectively completed, allowing us to draw 



1 Pamela is included in the Starlink distribution 'Hawaiki' and 
later releases. The Starlink Software Group homepage can be 
found at http://starlink.jach.hawaii.edu/starlink 
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Figure 2. Spectroscopic completeness of our sample. We show 
objects meeting our selection criteria for which there are Sloan 
spectra (black), targets for which we have taken our own spectra 
(grey), and targets yet to be observed (empty). Numbers above 
bins indicate the number of currently known emission line AM 
CVn binaries in our spectroscopy or the SDSS DR9 spectroscopic 
database. 



some preliminary conclusions. We estimate the number of 
systems in our sample, from the numbers we have found so 
far, by multiplying the numbers of AM CVn binaries found 
in each magnitude bin in Fig. [2]with the ratio of the total to 
observed number per bin. This suggests that there should be 
at least 5 more systems yet to be discovered in our sample, 
however, we note that with such small numbers there is a 
large uncertainty. 



3 CROSS-MATCHING 
3.1 GALEX 

The GALEX satellite has conducted a series of imaging and 
spec troscopic surveys, i ncluding the first UV all sky sur- 
vey (|Martin et al.l 12005). This consists of imaging in two 
broad bands, near- ultraviolet (NUV, 1770-2730 A) and far- 
ultraviolet (FUV, 1350-1780 A), of ~26,000 square degrees 
of sky, to a depth of 20.5 mag ( GALEX uses the AB photo- 
metric system). This is complemented by medium and deep 
imaging surveys with greater depths, covering smaller ar- 
eas. All candidates from the SDSS DR7 colour selection were 
matched against the GALEX Data Release 6 catalogue, tak- 
ing the closest neighbour within 2 arcsec. This matching ra- 
dius was chosen to ensure good coverage of the expected 
distribution of offsets between GALEX and SDSS astrom- 
etry, without intr oducing a significant number of false or 
multiple matches (|Morrissev et al.ll2007i ). 

We do not us e the predefine d GAL EX-SDSS cross- 
match described bv lBudavari et ab (2009), as we prefer the 
smaller FUV-NUV match radius used for t he standard cat- 
alogs , in order to minimize false matches jMorrissev et al.l 
I2007T I. We also prefer to have greater freedom in the choice of 
GALEX-SDSS matching radius. We expect that the results 
should be identical in most cases. 

We estimate the false match probability by repeating 
the search of the GALEX database after applying a ran- 
dom 10 arcsec offset to the coordinates of each SDSS source 
for which a GALEX counterpart has been found. The num- 
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Table 1. Summary of our observing campaign and the instrument setups used. The wavelength range and resolution achieved with 
each instrument are listed, as well as the number of candidates observed during each programme. Several of our candidates from 
SDSS DR7 have been observed by the SDSS as part of more recent data releases, those we have not observed previously ourselves 
are also listed here. 



Dates 


Telescope /Instrument 


Number of candidates observed 


Wavelength range (A) 


R at 5875A 


2008 Jan 09 - 2009 Apr 30 


Tillinghast/FAST 


225 




3500 


- 7400 


810 


2008 Feb 25 - 2008 Mar 02 


INT/IDS 


117 




3500 


- 8500 


630 


2008 Feb 27 


Keck-I/LRIS 


6 


3400 - 


5700, 6750 


- 9000 


730 


2008 June 06 - 2008 Sep 20 


VLT/FORS1 


18 




3600 


- 8600 


340 


2008 June 03 - 2009 Apr 25 


Hale/DBSP 


1 13 




3400 


- 8000 


540 


2008 Dec 25 - 2008 Dec 28 


WHT/ISIS 


68 




3200 


- 8200 


1160 


2009 Feb 26 - 2009 Apr 22 


Gemini-South/GMOS 


14 




3900 


- 6700 


640 


2009 July 21 - 2009 Aug 24 


Gemini-North/GMOS 


15 




3900 


- 6700 


640 


2009 Mar 17 - 2009 Mar 23 


INT/IDS 


95 




3500 


- 8500 


630 


2009 May 24 - 2009 May 28 


NOT/ALFOSC 


131 




3800 


- 9000 


180 


2009 June 18 - 2009 June 24 


WHT/ACAM 


205 




3800 


- 9200 


300 


2009 Oct 09 - 2009 Oct 15 


NTT/EFOSC 


59 




3800 


- 8000 


340 


2009 Nov 08 - 2009 Nov 14 


WHT/ACAM 


130 




3800 


- 9200 


300 


2010 Aug 19 - 2010 Aug 23 


WHT/ACAM 


88 




3800 


- 9200 


300 


2010 Nov 01 - 2010 Nov 06 


WHT/ACAM 


5 




3800 


- 9200 


300 


2011 Feb 02 - 2011 Feb 06 


NOT/ALFOSC 


55 




4000 


- 9000 


180 


2012 May 22 - 2012 May 23 


SOAR/Goodman 


5 




3800 


- 7000 


860 


2012 June 13 


Keck-I/LRIS 


1 




3100 - 


10000 


2100 


2012 July 13 


WHT/ISIS 


•1 


3800 - 


5200, 5600 


- 7100 


2200 


2008 July - 2009 June 


SDSS DR8 


17 




3800 


- 9200 


1800 


2009 June - 2011 July 


SDSS DR9 


11 




3650 - 


10400 


1800 



ber of matched coordinates divided by the total number of 
coordinates, gives an estimate of the false match rate; we 
derive a value of approximately 3.5 per cent. 

A total of 1622 of the 1947 SDSS DR7 objects have 
at least one measured UV magnitude, 1590 of these have 
a NUV detection, and 1138 are detected in both GALEX 
bands. 



3.2 UKIDSS 

The SDSS DR7 sample was also matched to the UKIRT 
(United Kingdom Infrared Telescope) Infrared Deep Sky 
Survey (UKIDSS) DR9 catalogue. When complete, the 
UKIDSS Large Area Survey will cover 4028 square degrees 
of sky, with imaging in four broad-band fil ters, Y, J, H an d 
K, to a F-band limiting magnitude of 20.3 (|Dve et alfc oOG). 
A matching radius of 1.0 arcsec was used for the UKIDSS 
catalogue (which has a greater astrometric accuracy than 
GALEX), this covers the majority of the offset distribution. 

Again we estimate the false match rate by offsetting 
the coordinates of the matched sample and repeating the 
search. We find a false match probability of ~1 per cent 
for the UKIDSS. Of our 1947 SDSS targets, only 516 are in 
the area that has been covered by UKIDSS, and 398 have a 
detection in at least one filter. 



4 SPECTROSCOPIC IDENTIFICATION 

We assign a spectroscopic classification to each object based 
on visual inspection of its spectrum. Tabl e [2] gives the nu m- 
bers of Cataclysmic Variables (CVs; see IWarnerlll995al for 
a review), A M CVn binaries, white dwarfs, subdwarfs (see 
lHeberll2009l for a recent review) and other objects identified 
in our sample. The white dwarfs are divided into several 



Table 2. Numbers of main object types identified in our sample. 
Objects classified as DA, DB, DQ and DZ white dwarfs have 
spectra dominated by hydrogen, helium, carbon and metal lines 
respectively; DC white dwarfs are those that exhibit a continuum 
spectrum. The subdwarf classification includes objects identified 
as sdB, sdOB or sdO. 



Class 


Number 


AM CVn 


7 


CV 


29 


quasar 


109 


galaxy 


2 


white dwarfs: 




DA 


120 


DB 


427 


DQ 


23 


DZ 


21 


DC 


73 


WD+dM 


1 


subdwarf 


184 


unknown 


407 


Total 


1403 



subclasses, DA white dwarfs are those with hydrogen atmo- 
spheres, DB those with helium atmospheres, DQ those with 
carbon dominated spectra, DZ those with spectra dominated 
by metal lines; and DC th ose with a continuum spectrum 
(see IWesemael et~all Il993l for an overview of white dwarf 
spectra). Fig.0shows example spectra of these main classes, 
with the features that identify them labelled. The AM CVn 
binaries, CVs and quasars are easily identified by their emis- 
sion lines; the systems that show only absorption in their 
spectra often have less certain classifications. A full list of 
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Figure 3. Identification spectra of select objects from our sample, demonstrating the appearance of the main object classes. The features 
that identify each object as a member of its class are labelled. The emission lines of helium and hydrogen present in the spectra of AM 
CVn binaries and CVs produce a strong contrast between them and the white dwarf contaminants, even at low S/N. The extremely 
broad, redshiftcd emission lines seen in quasar spectra, distinguish them from the other object types. 



candidates with classifications, coordinates, FUV, NUV and 
u, g, r, i, z magnitudes is given in Table IaTFI . 

Fig. [4] shows the He I 5875 against Ha equivalent width 
distribution of the spectra of our sample (here emission lines 
correctly have negative equivalent widths, this is ignored 
for simplicity elsewhere in the text). These were calculated 
using a Gaussian fit to the normalised spectrum, with fit 
constraints chosen to be consistent with the line fitting al- 
gorithm used by the SDSS. ff the fit failed, the offset from 
the expected central wavelength was larger than 20 A or the 
Gaussian dispersion was measured to be less than 0.5 A or 
greater than 100 A, the result was rejected and a value of 
zero taken. Objects appear in the expected regions of the 
diagram, with some scatter caused by noise. 

The AM CVn binaries (stars) lie close to the line 
EW(Ha) = A, with negative He I 5875 equivalent widths. 
SDSS J0804+1616 has significant He n 6559 emission that 
causes the apparent large equivalent width for Ha. The 
CVs (orange inverted triangles) are all found below the 

2 The full spectroscopic catalogue is available in electronic form 
at the CDS (http://cdsweb.u-strasbg.fr/). 



line EW(Ha) = EW(5875), and (with the exception of 
SDSS J203311. 78+134954.1, probably due to interstellar 
absorption) to the left of EW(5875) = A. As would 
be expected, the CVs identified all have hydrogen emis- 
sion, which is always stronger than helium emission. SDSS 
J121534.77+025726.6, identified as a white dwarf + M dwarf 
binary, is also found in this part of the diagram due to the 
hydrogen emission in its spectrum (most likely due to ir- 
radiation or chromospheric activity). Several quasars have 
apparently valid equivalent widths for some lines, and are 
scattered over the diagram; these are easily distinguished by 
visual inspection of their spectra, and also by their colour, 
which could be useful to distinguish those with low signal- 
to-noise spectra falling in the AM CVn region of Fig. 2] That 
there are no significant large deviations from the expected 
position for any object type, suggests that our classification 
by visual inspection was successful. 

There are 407 spectra which have no clear identifiable 
features, so no likely classification has been assigned. These 
objects are generally fainter, resulting in lower signal-to- 
noise ratio spectra, which makes weak features more diffi- 
cult to discern. The majority of these objects are likely white 
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Figure 4. Equivalent width of He I 5875 vs Ho for our sample, the right hand panel shows the crowded area of the diagram, occupied 
mainly by white dwarfs, in more detail. The black stars and orange inverted triangles represent AM CVn binaries and CVs. Quasars, 
DB white dwarfs, other types of white dwarf, subdwarfs, and WD + M dwarf systems are represented by red triangles, cyan squares, 
blue diamonds, green circles and black squares respectively. Grey dots indicate those candidates that could not be classified. Equivalent 
widths have been calculated using a Gaussian plus sloped continuum fit; the large values for some unclassified objects are largely the 
result of noise. Objects appear in the expected regions with allowance for some scatter caused by noise. Note that He II 6559 emission 
would not be distinguished from Ha by the fitting algorithm. Several quasars have apparently valid equivalent widths for some lines, and 
appear randomly distributed, but are distinguishable using other spectral features. 



dwarfs with weak lines or DC white dwarfs; we do not assign 
this classification as we do not want our assumptions as to 
their nature to affect our understanding of the colour space. 
With the exception of 18 candidates, none of these objects 
have helium emission with an equivalent width greater than 
3 A as the dominant feature in their spectra, ruling them 
out as possible AM CVn binaries. The measured equivalent 
widths for several of these exceptions that have possible he- 
lium emission, have very large errors, and the measurements 
are rejected on this basis. Fourteen of these unclassified ob- 
jects are not consistent with having an He I equivalent width 
of A; although the apparent emission would be classified as 
noise by visual inspection, they cannot be ruled out as AM 
CVn binaries and should be re-observed to obtain higher 
signal-to-noise spectra. 

The AM CVn binaries discovered in the SDSS spec- 
troscopic database, and via this survey, all have equivalent 
width > 13 A for one of the 4686 A or 5875 A helium 
emission lines. This indicates our ability to identify can- 
didate AM CVn binaries matches th at of previous searches 
of the SDSS spectroscopic database llAnderson et all 120051 ; 
iRoelofs et al. 1 120051 ; lAnderson et al.l l2008t ); and we can be 
confident that we would be able to detect the AM CVn bin a- 
ries expected from predictions (jRoelofs et al.|[2007d .[2009). 

It must also be noted that high-state AM CVn bina- 
ries - those with short orbital periods and those in out- 
burst - typically show absorption, rather than emission, 
in their spectra, due to the optically thick accretion disc. 
The spectra of these systems look very similar to those 
of DB white d warfs, AM CVn itself w as originally iden- 
tified as such (Greenstein fe Matthews! 119571 ). there are 



however, some differences in the line strengths that allow 
them to be distinguished with high-quality s pectra (e.g. 
iRobinson fc FaulkneJl975l ; lFontame et al . 2011). At the low 
spectral resolutions and signal-to-noise ratios employed in 
this survey, these high-state systems are essentially indis- 
tinguishable frojii_rjB_wdTitedwarfs, or possibly DC white 
dwarfs (e.g. lO'Donoghue et al.lll987l ). 

The short period AM CVn binaries are expected to be 
far less numerous at the high Galactic latitudes of the SDSS, 
as the evolution from t urn-on of mass transfe r to P OI b > 30 
min is relatively rapid IjNelemans et al.ll200lr ). Using a sim- 
ple model for the evolu tion of a white dwarf channel system 
(jNelemans et al.ll200il ). we find that an AM CVn binary 
with t ypical values for th e mass at period minimum (fol- 
lowing [Levita^eF|al][20l3), has an orbital period below 30 
minutes for 5.8 per cent of its life between period minimum 
and P or b = 60 minutes. The high state accretion discs in the 
short period systems also make them much brighter, result- 
ing in increased numbers in a magnitude limited sample. 
With an estimate of both the relative numbers and their 
absolute magnitude as a function of period, it is possible 
to calibrate the AM CVn space density b ased only on the 
emission line systems jRoelofs et al.ll2007cf ). 

4.1 Cataclysmic variables 

Fig. [5] shows the spectra of the CVs identified in our sample 
via their hydrogen emission. The range of line strengths and 
ratios seen in Fig. [4] are clear. The broad double-peaked line 
profile characteristic of an accretion disc is evident in several 
objects. 
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Figure 5. Identification spectra of the CVs discovered in our sample, displayed with 2 pixel Gaussian smoothing. Dotted lines indicate 
the wavelengths of the strongest hydrogen lines. 
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Figure 6. Equivalent width of He I 5875 vs Ho for our CV sample, 
and the Sloan CV population. Symbols have the same meanings 
as in Fig. [4] and the crosses represent the Sloan CVs. 



SDSS J074859.54+312512.7, SDSS J075107.51+300628.5 
and SDSS J131432. 11+444138 .8, were previously identified 
as dwarf nova candidates by IWils et all |20ic!) based on 
their variability. The spectra presented here confirm these 
systems as CVs. SDSS J225417.54+074227.3 (USNO- 
A2.0 0975-21 112378) was p r eviou sly identified as a QSO 
candidate bv lAtlee fc Gouldl (|2007i ). 

We compute the He n 4686/H/3 line ratios for the 
CVs with significant He II emission, in order to iden- 
tify possible magnetic or nova-like systems. The five 
systems with an equivalent width ratio greater than 
0.5 are listed in Table [3] The He n 4686 equivalent 
widths measured for SDSS J175320.59+251649.1 and SDSS 
J212617.62+192320.1 may be affected by blending with 
nearby lines, but are still large compared with most CVs. 
The weak lines exhibited by SDSS J003719.29-215714.4 and 
SDSS J175320.59+251649.1, combined with the strong He n 
suggest they may be high accretion rate nova-like systems. 
SDSS J111126. 83+571238.6 also shows unusually strong He I 
emission and is discussed further below. 

Fig-E]shows tjjg equivalent width distribut ion of our CV 
samp le from Fig. [4J with the Sloan CV sample (|Szkodv et al.l 
l201ll l added for comparison. The distributions of the two CV 
populations look similar. The strongest Ha emitters fall out- 
side our colour box, and none of the CVs in our sample have 
Ha emission as strong as the strongest Sloan CVs. There 
are two obvious outliers, SDSS J111126.83+571238.6 (SBSS 
1108+574) from our sample a nd SDSS J112253.3 -111037.6 
(CSS100603:112253-111037, IBreedt et all |2012| ) from the 
Sloan population, that have much stronger helium emission 
relative to hydrogen than the majority of CVs (this is also 
clear from the spectrum of SDSS Jllll+5712 shown in Fig. 
[5J. These may represent hybrid CV - AM CVn binaries, 
and may be AM C Vn binaries forming via the evolved CV 
form ation channel ( Podsiadl owski et al.ll200"3l ; IBreedt et al.l 
120121 ). 



4.2 DQ white dwarfs 

Fig.[7]shows the spectra of the DQ white dwarfs identified in 
our sample. These objects have been identified as DQs due 
to the presence of either neutral carbon lines or molecular 
C2 Swan bands in their spectra. A number also show weak 
Ha absorption and are classified as DQA. 

The expected correlation of the transition from atomic 
to mole cular carbon f eature s with colour (or temperature), 
seen by lHarris et"aH |2003l ) in SDSS DQ white dwarfs, is 
also seen in our DQ sample. 



4.3 DZ white dwarfs 

Fig.[8]shows the spectra of the DZ white dwarfs identified in 
our sample. These objects have been identified as DZs due 
to the presence of the Ca 11 H and K lines in their spectra. 
A number also show weak Ha absorption and are classified 
as DZA. 

SDSS J105338.16+285245.6 (USNO-A2.0 1125- 
06201629) was previo usly identified as a QSO candidate by 
lAtlee fe Gouldl (|2007l 1. 

All the DZs in our sample show calcium absorption 
with no detection of other meta ls, similarly to the major- 
ity of SDSS DZ white dwarfs (|Harris et all 120031 ). They 
are all cool in (u — g, g — r), with the exception of SDSS 
J132430. 43+055316.0, which also shows significant hydro- 
gen in its spectrum. 



5 RESULTS 

5.1 Identification spectra of two new AM CVn 
binaries 

The extracted spectra of SDSS J104325. 08+563258.1 (here- 
after SDSS J1043) and SDSS J173047.59+554518.5 (here- 
after SDSS J1730) are shown in Fig. [5] The helium emission 
lines at 4471, 5875, 6678 and 7065 A are clearly detected, 
and no hydrogen is present. This suggests that SDSS J1043 
and SDSS J1730 are likely new members of the AM CVn 
class, although further observations are required to confirm 
their ultra-compact binary nature and determine their or- 
bital periods. 

SDSS J1043 was identified as a candidate dwarf nova in 
th e second GALEX Ultr aviolet Variability catalog (GUVV- 
2; IWheatlev et al.l [2008). This catalog consists of objects 
identified as varying by more than 0.6 mag in multiple NUV 
observations of ~161 deg 2 of the sky. The repeat observa- 
tions of SDSS J1043 show two outbursts of >2 mag dur- 
ing the two year win dow. The Catalina Real-Time Tran- 
sient Survey (CRTS; iDrake et all 120091 ) also shows three 
outbursts of >2 mag during a ~4.5 year time span. 

The quiescent spectrum of SDSS J1043 shows the char- 
acteristic helium emission of AM CVn binaries, but is un- 
usual in that it shows absorption wings around these lines. 
Simil ar spectra are shown by V406 Hya (P or b = 33.8 min- 
utes; IRoelofs et all l2006al) a nd SDSS J1240 (P or b = 37.4 
minutes; IRoelofs et al.l 120051 ) . in which the DB absorption 
featu res are assumed to be due to the accretor (|Roelofs et al.1 
120041 ; which is still hot compared to the longer period AM 
CVn binaries that show no absorption features). 

TableUlists the equivalent widths of the main lines seen 
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Table 3. He II 4686 and H/3 emission line equivalent widths and line ratios for CVs 
with strong He II emission. 



Name 




EW He II 4686 (A) 


EW H/3 (A) 


Line ratio 


SDSS J003719. 29-215714.4 
SDSS J111126.83+571238.6 
SDSS J175320.59+251649.1 
SDSS J204643.30-000630.2 
SDSS J212617.62+192320.1 


-3.3 ± 0.1 
-19.0 ± 2.0 
-6.0 ± 0.5 
-3.8 ± 0.3 
-18.0 ± 2.0 


-3.8 ± 0.1 
-34.0 ± 2.0 
-6.9 ± 0.5 
-5.9 ± 0.3 
-33.0 ± 2.0 


0.87 ± 0.03 
0.56 ± 0.07 
0.87 ± 0.10 
0.64 ± 0.06 
0.55 ± 0.07 
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Figure 7. Identification spectra of the DQ white dwarfs discovered in our sample, displayed with 2 pixel Gaussian smoothing. Dotted 
lines indicate the wavelengths of the strongest lines of neutral carbon, dashed lines indicate the C2 Swan band heads. The dot-dashed 
line indicates the wavelength of Ha, which is present in several objects. SDSS J1120— 1102 was observed as part of SDSS DR8, and has 
not been observed separately by us. 
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Figure 8. Identification spectra of the DZ white dwarfs discovered in our sample, displayed with 2 pixel Gaussian smoothing. Dotted 
lines indicate the wavelengths of the calcium H and K lines. The dot-dashed line indicates the wavelength of Ha, which is seen in several 
objects. 



Table 4. Equivalent width and Gaussian F WHM for the prominent emission 
lines in SDSS J1043 and SDSS J1730. Estimated errors are largely due to 
uncertainty in the continuum. 



Line 



SDSS J1043 
EW (A) FWHM (km s" 1 ) 



SDSS J1730 
EW (A) FWHM (km s" 1 ) 



He I 4387 






-9 ± 1 


1100 ± 100 


He I 4471 


-4 ± 1 


3300 ± 400 


-20 ± 1 


1100 ± 100 


He II 4686 


-12 ± 2 


5000 ± 500 


-21 ± 1 


2600 ± 200 


+ He I 4713 










He I 4921 






-13 ± 1 


820 ± 20 


He I 5015 


-10 ± 2 


3200 ± 500 


-19 ± 1 


1010 ± 30 


He I 5875 


-16 ± 2 


2900 ± 200 


-53 ± 2 


820 ± 10 


He I 6678 


-10 ± 2 


3000 ± 300 


-45 ± 2 


730 ± 10 


He I 7065 


-20 ± 4 


3500 ± 300 


-40 ± 2 


630 ± 40 



Values marked '— ' could not be measured reliably. 
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Figure 9. Identification spectra of SDSS J1043 (taken with LRIS on Keck I on 2012 June 13) and SDSS J1730 (taken with ACAM on 
the WHT on 2009 November 14). The wavelengths of the strongest He I and hydrogen lines are indicated by the dotted and dashed lines. 
The lack of hydrogen features in the spectra is clear. SDSS J1043 shows clear double-peaked line profiles, in contrast to the narrow lines 
displayed by SDSS J1730, this is likely due to a lower inclination of SDSS J1730 . SDSS J1043 also shows absorption wings around the 
He I 4026 and 4471 A lines, thought to be due to the accretor (e.g. iRoelofs et alj|2004tl . 



in the spectra of SDSS J1043 and SDSS J1730. The strong 
helium lines in SDSS J1730 are com parable to tho se shown 
by SDSS J1411 (P orb = 46 minutes; lRoelofsll2007l) and GP 
Com (P orb = 46.5 minutes; iNather et al.lll98ll ). Together 
with the absence of helium absorption at shorter wave- 
lengths, this suggests that SDSS J1730 may be at the long 
period e nd of the AM CV n period distribution (P or b >~ 40 
minutes. iRau et al.l 120101 ). Whilst the equivalent widths of 
the lines are similar to those shown by other known AM CVn 
binaries, the lines are narrow compared with those shown by 
other systems in our sample. The small FWHM of the emis- 
sion lines in SDSS J1730 may indicate a low inclination disc. 
The lack of obvious double peaked structur e to the lines also 
sugge sts a low inclination for this system |Horne fc Marshl 
Il986h . 



5.2 AM CVn equivalent width — period relation 

The orbital periods and He I 5875 equivalent widths for AM 
CVn binaries with emission line spectra are collected in Ta- 
ble [5] objects for which no spectrum was available and no 
measured equivalent width could be found are excluded. Fig. 
1 101 shows the correlation of the equivalent width of the He I 
5875 emission line with orbital period for AM CVn binaries. 
There is a general trend for larger equivalent widths with 
longer binary periods. This is probably due to the drop in 
continuum flux, against which the lines are measured, due to 
the cooling of the accreting white dwarf, and drop in accre- 
tion continuum as the mass transfer rate decreases towards 
longer periods. The expected increase in disc size at longer 
periods, as the orbit expands, may also contribute to the 
effect. 

This trend would not be expected to hold below a period 
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Figure 10. Equivalent width of the He I 5875 line versus orbital 
period for AM CVn binaries. The black circles represent those AM 
CVn binaries discovered from the SDSS. The dotted lines indicate 
He I 5875 EWs for the AM CVn binaries discovered via this survey 
whose periods are unknown. The dashed line and shaded area 
indicate our detection threshold. 



of ~20 minutes, where the disc is expected to be in a stable 
high state. ES Cet is included as it shows helium in emission, 
although it might not be expected to follow the same trend 
due to its much shorter orbital period. ES Cet in fact shows 
very strong He II 4686 emission, but the weaker 5875 A line 
appears to match the trend of the longer period systems 
reasonably well. 

Whilst there is considerable scatter that prevents an 
accurate estimate of the orbital period based only on the 
equivalent width, the relation does suggest a period range 
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Table 5. Orbital period and equivalent width of the He I 5875 line for the AM CVn binaries with 
quiescent spectroscopy. 
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in which an object might be expected to be found. Based on 
Fig. 1101 SDSS J1043 is expected to have an orbital period 
between 25 and 45 minutes; SDSS J1721 is expected to have 
an orbital period in the range 25-55 min; and SDSS J1730 is 
predicted to have an orbital period between 40 and 55 min. 

This trend also demonstrates that AM CVn binaries 
with orbital periods in excess of 30 min should have emis- 
sion lines sufficiently strong to be well above our detection 
threshold. We note that the drop in continuum flux at longer 
orbital periods is expected to make the longer period sys- 
tems fainter, and less likely to be found within our g-band 
magnitude cut, or the SDSS spectroscopy. This may intro- 
duce a selection bias, only the systems with very strong 
emission appear sufficiently bright to have been observed, 
and hence there is the possibility of an undetected popula- 
tion of long period, low equivalent width systems. 

5.3 Colours of the sample 

The (u — g, g — r) and (g — r, r — i) colour diagrams of 
the sample and the colours of a subset of the known AM 
CVn binaries (black stars) are shown in Fig. 1111 The most 
significant of the classifications from Table [2] are also indi- 
cated with different markers. The seven systems discovered 
in the observed part of the sample further constrain the re- 
gion of colour space occupied by AM CVn binaries. The 
colours of AM CVn itself may be affected by saturation due 
to its brightness, it has therefore been excluded, u — g val- 
ues calculat ed from the UBV p hotometry of CR Boo given 
in table 1 of IWopd et al] (1 19871 ) (using conversions given by 
IJordi et alj|2006h . suggest that the colours can vary signif- 
icantly over an outburst cycle. CR Boo was likely to have 
been in an intermediate state at the time of the SDSS pho- 
tometry (g — 15.6 compar ed to g ~ 17.1 i n quiescence, 
based on photometry from IWood et al.l 1 19871 ) , which may 
have caused the unusual u — g colour. 



The initial selection of the sample l|Roelofs et al.ll2009l ) 
was based on the colours of the nine emission line AM CVn 
binaries with SDSS photometry known at the time. Fig. [11] 
shows the colours of 24 of the 35 currently known AM CVn 
binaries. Whilst many of the systems discovered since have 
been found through this survey, all long period systems, and 
most intermediate period systems with quiescent SD SS pho- 
tome try, fall within our colour box (see also iLevitan et al.l 
|2012| . for a discussion of the colours of outbursting systems 
found by the Palomar Transient Factory). It should be noted 
that the survey targets the specific area of colour-colour 
space in which the AM CVn binaries in the SDSS spectro- 
scopic database are found. We are therefore only capable of 
finding systems that share these colours, and any long pe- 
riod AM CVn binaries that lie outside this colour box will 
remain undetected. The biases in the S DSS spectro scopic 
database are discussed bv lRoelofs et~aH (|2007d . l2Q09h . they 
conclude that it is unlikely that a large fraction of AM CVn 
binaries lie outside this region. This more robust knowledge 
of the location in colour space of AM CVn binaries, and 
other objects in this region, allows us to improve upon the 
initial selection criteria and reduce the sample size. 

It is clear from Fig.[TT]that the known AM CVn binaries 
with SDSS photometry occupy the region with 

g-r < -0.1; 

approximately 81 per cent of the quasars (red triangles) have 
g — r above this limit, and many could be safely removed. It 
should be noted that the spectroscopic completeness of the 
SDSS increases with increasing g — r over the densely pop- 
ulated area of this region, where objects have been targeted 
for quasars (Fig.[TJ. The original sample of systems from the 
SDSS spectroscopic database is thus slightly biased towards 
the red cut-off. This makes it unlikely that a large fraction 
of AM CVn binaries lie beyond g — r = — 0.1, and we con- 
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Figure 11. Colours of the known AM CVn binaries (black stars) together with our sample of candidates from SDSS DR7. The filled 
stars are the confirmed longer period systems, and the AM CVn binaries discovered via their emission line spectra in the SDSS or this 
survey. The open stars represent those systems to which we do not necessarily expect to be sensitive, from left to right (left hand panel) 
or from top to bottom (right hand panel), these are HM Cnc (right hand panel only), SDSS J1908, SDSS J2047, PTF1 J2219, ES Cet, 
PTF1 J0857, PTF1 J0943, PTF1 J0435 (left hand panel only) and CR Boo. The orange inverted triangles represent CVs. Quasars, DB 
white dwarfs, other types of white dwarf, subdwarfs, and WD + M dwarf systems are represented by red triangles, cyan squares, blue 
diamonds, green circles and black squares respectively. Grey dots indicate those candidates that could not be classified, and black dots 
those that have not yet been observed. The solid line shows the blackbody cooling track, the dotted and dot-dashed lines represent the 
DA and DB white dwarf cooling tracks, and arrows represent reddening vectors for an extinction A(g) = 0.2. The dashed line represents 
the new colour cut. 



sider the risk to our completeness to be small in comparison 
to the gain in efficiency. 

The addition of UV photometry allows further examina- 
tion of the sample in additional colour spaces. Fig. ll2l shows 
colour diagrams of the part of the sample with GALEX 
NUV detections, these diagrams include 17 of the known 
AM CVn binaries. UV e xtinction was estima ted using the 
Rv = 3.1 prescription of lCardelli et al. and the full 

Galactic reddening according to Schlegel et al. fl998), in or- 
der to match the procedure used for the SDSS. Note that 
CR Boo and SDSS J1043 were much closer to their max- 
imum brightness states when observed by GALEX, com- 
pared with their SDSS photometry, causing the large neg- 
ative NUV — u. There are also several other objects with 
NUV — u < —1, for which different brightness states are a 
possibility (NUV — u ~ — 1 is the limit for a hot blackbody). 

The NUV — u colour greatly increases the separation of 
some object types compared with only Sloan colours. The 
subdwarfs and white dwarfs occupy mostly separated areas 
of the diagrams in Fig. 1121 NUV — u also separates the subd- 
warfs and the quasars. Many of the subdwarfs (green circles) 
in the upper right corner of the left hand panel of Fig. [T^] 
could also be discarded, as the AM CVn binaries do not 
spread into this region. By requiring 

NUV-u > 4.34( 5 -r) +0.5, 

46 per cent of the subdwarfs can be removed, and the sample 
size is reduced by 11 per cent. 

The DB white dwarfs (cyan squares) lie around and 



close to the blackbody cooling track in (NUV — u, r — i) 
colour space, whereas the AM CVn binaries all lie above. 
With a cut close to this cooling track, and parallel to the 
reddening vector, 

NUV -u< 6.76(r + 1.85, 

approximately one fifth of the DB white dwarfs can be re- 
moved from the sample, significantly reducing the numbers 
of the main contaminant, and reducing the total size of the 
sample to ~1500. 

We emphasize that these additional cuts increase the 
risk of missing some AM CVn binaries that are present in 
the SDSS photometry, making it more difficult to judge our 
completeness, and the accuracy of our space density. How- 
ever, the increase in efficiency offers a significant advantage 
in reducing the observing time required to complete such a 
survey. 

The UKIDSS colours of the sample are shown in Fig. 1131 
The separation of quasars from the rest of th e sample is clear 
even with the limited sky coverage (see also iRichards et al.l 
l2009l ; IWu fcjiall2010l ; for a discussion of quasar selection via 
infrared colours). As only four of the known AM CVn bina- 
ries have been covered by UKIDSS, and none with sufficient 
detections whilst in the same brightness state to appear in 
Fig. 1131 any attempt to remove targets from the sample 
based on infrared colours would be unreliable. Greater cov- 
erage with infrared photometry could, however, allow fur- 
ther improvements to the sample size. The WISE mission 
( Wrig ht" et al.l l2010) has recently conducted an all sky in- 
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Figure 12. Colours of the known AM CVn binaries (black stars) together with our sample of candidates from SDSS DR7 combining 
GALEX NUV detections. The symbols have the same meanings as in Fig. 1111 The open stars represent, from top to bottom, CR Boo, 
PTF1 J2219, ES Cet, HM Cnc (right hand panel only) and PTF1 J0857. The large negative NUV — u colour of CR Boo and SDSS J1043 
is caused by their varying brightness states between the epochs of GALEX and SDSS observations. The solid line shows the blackbody 
cooling track, the dotted and dot-dashed lines represent the DA and DB white dwarf cooling tracks, and arrows represent reddening 
vectors for an extinction A(g) = 0.2. A significant portion of DB white dwarfs lie below of slightly above the blackbody cooling track in 
(NUV — u, r — i) colour space, whereas the AM CVn binaries all lie above this track. The dashed lines represent the new colour cuts. 
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Figure 13. Infrared colours of our sample. The symbols have the same meanings as in Fig. 1111 Both diagrams show the potential for 
removing quasars that infrared data provides, however, greater coverage is required to enable reliable cuts to be made. 



frared survey, however, the depth of the survey is insufficient 
for the majority of known AM CVn binaries. 



5.4 The AM CVn binary space density 

The discovery of only seven AM CVn binaries in the 72 
per cent of the sample that has been observed, taking into 
account our spectroscopic completeness (Fig. indicates 
the presence of significantly fewer AM CVn binaries in 
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the SPSS photometr i c data base than initially estimated by 
IRoelofs et all i|2007d . 120091 ) (~40 new AM CVn binaries in 
the original sample of ~1500 candidates from SDSS DR6; 
there are ~2000 objects in the DR7 sample). This suggests 
that the space density of AM CVn binaries is lower than 
previously predicted, or that there is some problem with our 
understanding of the biases in our sample or t he population 
models . Base d on the analysis presented in IRoelofs et al.l 
l|2007d . [200gh . we consider the former to be likely. 

Fig. [2] shows that our follow-up is almost complete to a 
depth of g — 19. We use this essentially complete sample of 
the brighter objects in the AM CVn colour box to adjust ear- 
lier space density estimates, although the large uncertainty 
resulting from such small samples must b e noted. Following 
the prescription of IRoelofs et al.l (|2007d ). we calculate the 
expected magnitude distribution of the AM CVn binaries 
in the SDSS photometric database. The number of systems 
expected is then compared to the number found in the SDSS 

DR7 area. 

B ased o n the optimistic model from iNelemans et al.l 

(|200ll . 12003 ). and the six AM CVn binaries discovered in 
the SDSS spectroscopy, IRoelofs et al.l (|2007d ) calculated an 
observed space density 1.5 x 10 -6 pc -3 . This corresponds to 
an expected 35 AM CVn binaries to a depth of g = 20.5 (see 
Table [6]), in the DR5 photometric database (note that this 
number is ess entially independent of the formation channel; 
see table 1 of IRoelofs et al.ll2007d ). To find the distribution 
of systems the p opulation synt h esis pr edicts for our sample, 
we multiply the IRoelofs et al.l (|2007d ) distribution by the 
ratio of the photometric area of DR7 (11663 deg 2 ) to that 
of DR5 (8000 deg 2 ), ~1.46. This gives expected numbers of 
11 systems with g ^ 19, and 51 systems with g ^ 20.5. 

We use the 4 long period emission line AM CVn binaries 
found in our essentially complete g ^ 19 sample to scale 
the expected numbers and observed space density (Tabled . 
This gives the total number of AM CVn binaries expected 
in the DR7 photometry (g ^ 20.5) as 18, corresponding to 
a space density for AM CVn binaries of 5 x 10~ 7 pc -3 . 

The small sample of 4 known AM CVn binaries with 
g ^ 19 contributes an intri nsic uncertainty of 50 per cent 
to the derived space density. IRoelofs et ail l|2007d ) estimate 
the uncertainties in the parametrization of temperature and 
absolute magnitude with orbital period as leading to ~10 
per cent and ~32 per cent variations in the result. Combin- 
ing these contributions leads to an estimated 60 per cent 
uncertainty in our value for the space density. 

The space density calculations assume models for the 
intrinsic magnitude distribution of the AM CVn binaries 
(|Roelofs et al.ll2007d ). It is clear from Fig.[2]that we find an 
increasing number of these systems towards fainter magni- 
tudes. This may reduce the discrepancy between the previ- 
ously predicted space density and the numbers of systems 
found so far in our sample, but it may also indicate a signif- 
icant deviation from the modelled magnitude distribution. 



6 DISCUSSION 

The results from our survey so far significantly improve our 
knowledge of the colours of AM CVn binaries, and the other 
objects found in the same region of colour space. The seven 
new AM CVn binaries lie in the central area of the selected 



colour region, implying that the initial selection should in- 
deed have contained the majority of AM CVn binaries, but 
also that it could have been more efficient. The reddest part 
of the g — r distribution is largely made up of quasars and 
unclassified objects, and can be removed from the sample. 

Cross matching the sample with the GALEX catalogue 
provides a further colour space that proves to be very useful 
in separating different classes of object. The majority of the 
subdwarfs in our sample are brighter in NUV relative to 
u than the white dwarfs, causing the gap between the two 
populations that is not seen with only Sloan colours. That 
the AM CVn binaries are slightly redder than many of the 
DB white dwarfs in r — i, perhaps due to the contribution 
of the accretion disc, is very useful when combined with 
their distribution in NUV — u. This increased splitting of 
DB white dwarfs from the other objects in the sample in 
(NUV — u, r — i) colour space, allows about one fifth of 
this major contaminant to be safely discarded. Subdwarfs 
and AM CVn binaries become increasingly separated with 
increasing g—r in (NUV— u, g—r) colour space, such that the 
majority of the redder subdwarfs can be cut. Making these 
cuts parallel to the reddening vectors ensures that objects 
are not removed unintentionally due to poorly estimated 
extinction. The small number of objects with NUV — it < — 1 
or NUV — it > 1.5 are kept, as the large values may be caused 
by varying brightness state between the epochs of GALEX 
and SDSS observations. 

Combining the two new cuts involving NUV — u, and 
the more efficient g — r cut-off, in addition to the original 
u ~ 9i 9 ~ r an d r — i colour cuts used to produce the sam- 
ple, reduces the total size by 43 per cent. Removing 268 of 
the 544 remaining targets, brings the goal of detecting all of 
the hidden AM CVn population much closer to completion. 
These extra cuts, whilst greatly increasing our efficiency, 
risk reducing our completeness. Although the new colour 
criteria avoid the regions in which AM CVn binaries have 
been found, it is possible that there exists some as yet unde- 
tected part of the population that deviates from the current 
distribution. Any such systems would not be expected to 
represent the majority of the population, and we consider 
the risk in losing them to be small. Increasing the efficiency 
of our survey becomes increasingly important as we move to 
fainter targets, where greater numbers of AM CVn binaries 
have been found (see Fig. [2}, in order to keep the project 
feasible. 

Fig. H4l shows the spectroscopic completeness of the sur- 
vey colour box, combining both SDSS spectroscopy and our 
own. Our spectroscopic followup has been conducted with- 
out any specific colour bias beyond the initial selection cri- 
teria. However, the magnitude bias shown in Fig. pleads to 
a slight bias towards the bluer part of the (u — g, g — r) 
distribution. The lower density of objects at the edges of 
our survey area leads to a significantly higher completeness. 
There is no clear bias that would affect the AM CVn distri- 
bution we observe. 

The observed trend of increasing equivalent width of 
emission lines with orbital period, coupled with the expec- 
tation that AM CVn binaries will evolve rapidly to longer 
periods, allows us to be confident that we are able to de- 
tect the majority of the population. This is important as we 
intend to use the number we find to estimate the space den- 
sity of the entire population. The equivalent widths of the 
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Table 6. Expected numbers of AM CVn binaries in the AM CVn colour box and corresponding space 
densities. Note that iRoelofs et al]|2007d use a depth of g = 21, whereas our survey has a limit of 
g = 20.5. 



Distribution 


Expected number of AM CVn 
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Space density (pc 3 ) 
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21 
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53 
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Roelofs et al. (2007c) DR7 area 


11 51 


77 
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Figure 14. Completeness of the survey as a function of colour, 
u — g and g — r, to a limiting magnitude g = 20.5 (dereddened), 
including both SDSS spectroscopy and our own. The AM CVn 
binaries with spectra taken as part of the SDSS or this survey, 
are indicated by star symbols. The solid line marks the blackbody 
cooling track, the dotted and dot-dashed lines indicate model 
cooling sequences for DA and DB white dwarfs. 

SDSS AM CVn binaries with unknown periods suggest that 
many of them are at the shorter period end of the di stribu- 
tion that we would expect (see IRoelofs et al.l l2007d ) . This 
suggests that there should be more systems with P or b > 50 
min yet to be discovered in our sam ple, as this is where 
AM CVn binaries sh ould accumulate (jNelemans et al.ll200ll ; 
iNissanke et~aT1l2012l '). 

Using our essentially complete g ^ 19 sample, we cal- 
culate the expected number of AM CVn binaries in SDSS 
DR7 to be 18 ± 9 {g «S 20.5). Since we know of 15 longer 
period systems with SDSS photometrjjf], this suggests that 
there should be 3 undiscovered systems left in our sample, 
though the considerable uncertainty on this number leaves 
us unable to make firm predictions. From the spectroscopic 
completeness as a function of g-band magnitude (Fig. [2}, 
we would expect to find ~5 more systems in the remainder 
of our sample, this is consistent with the number derived 
from our revised space density estimate, but also has a large 
uncertainty. 

3 This number includes the 6 systems discovered in the SDSS 
spectroscopy, the 7 systems discovered via our survey, CP Com, 
and SDSS J1427-0123 - which has not been observed spectroscop- 
ically by us or the SDSS. 



IRoelofs et all (|2007d ) calculated the AM CVn space 
density to be 1.5 x 10~ 6 pc -3 , based on the six systems 
found in the SDSS spectroscopic database. Our revised value 
is 3 times lower than thi s, and 50 times l ower t han the op- 
timistic model from the iNelemans et al.l (|200ll ) population 
synthesis predicts. 

INissanke et all (|2012T l suggest that the lower space den- 
sity could be explained by a chang e to the galactic disc 
model used by INelemans et all (|200ll ). The distances of AM 
CVn binaries from our solar system, and hence the ampli- 
tude of their gravitational wave signal, depends on the as- 
sumed scale height for AM CVn binaries in t he Galaxy. In 
the proposed alteration (Nissankc et al.l [20121 '). the old sys- 
tems, that are the ones we would detect as emission line 
systems in the SDSS sample, have a larger scaleheight, and 
so are farther away. However, the young systems, that are 
the majority of those that would be detected via their grav- 
itational wave e mission, remain unchanged from the origi- 
nal distribution (|Nelemans et alj|200ll . l2004t ). This results in 
the same total numbers of AM CVn binaries in the Galaxy, 
but the long period systems are at greater distances, and 
are hence fainter. This could expla in why optical surveys 
find less systems than pr edicted by INelemans et al.l l|200lf ) 
and IRoelofs et al.l (|2007d ). but has no significant effect on 
the numbers of AM CVn binaries detectable by space-based 
gravitational wave missi ons (compare cases 1 and 5 shown 
in table 2 and table 3 of INissanke et al.ll2012l ). 

Alternatively, the lack of systems could be explained 
as problems with the binary population synthesis; if fewer 
systems than expected survive to become stable mass trans- 
ferring AM CVn binaries, the predicted space density would 
fall. It should be noted that there is a large intrinsic un- 
certainty in the popu lation synthesis numbers, this is why 
INelemans et alj (|200ll ) have optimistic and pessimistic mod- 
els. It is also possible that the discrepancy is so large that 
it can only be explained by changes to both the popula- 
tion synthesis and the expected brightness distribution of 
the AM CVn binaries. 

This survey is based on the colours of the known emis- 
sion line AM CVn binaries, and hence targets the long pe- 
riod systems (P or b > 30 min). Whilst the majority of all 
systems fall within our colour box, we are only sensitive to 
intermediate period systems in their low state. In their high 
state these AM CVn binaries normally appear much like DB 
white dwarfs, and if they were also in their high state at the 
time the SDSS observed them, we would be unable to recog- 
nise them as AM CVn binaries. It should be noted that a 
small percentage of the total AM CVn population should be 
found as outbursting systems. We estimate that AM CVn 
binaries spend less than 6 per cent of their lifetime in this 
state. These frequent outbursting systems (20 < P or b <~ 30 
min) can be detected more efficiently from their variability, 
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using synoptic surveys, as in the case of PTF1 J0719+4858 
(|Levitan et alj|201ll . l20l3 >. Follow-up of these systems, how- 
ever, may be more challenging as they are only expected to 
be found in significant numbers by increasing the depth of 
the search. SDSS J1043 is particularly interesting as it repre- 
sents the overlap between these two methods, detected both 
as a result of its colours and its variability. These compli- 
mentary methods mostly sample different parts of the AM 
CVn orbital period distribution, and together, will lead to 
a better understanding of both this period distribution, and 
the space density. 



7 CONCLUSION 

We have presented an update on the status of our spec- 
troscopic survey aimed at uncovering the expected hidden 
population of AM CVn binaries in the SDSS photometric 
database. The results currently indicate a lower space den- 
sity than predictions suggest. Based on the brighter part 
of our sample, we calculate an observed space density of 
(5 ± 3) x 1CT 7 pc" 3 . 

We have reported the discovery of two candi- 
date AM CVn binaries found via this survey, SDSS 
J104325.08+563258.1 and SDSS J173047.59+554518.5. 
SDSS J1043 exhibits the helium absorption and low equiv- 
alent width emission lines shown by the AM CVn binaries 
with orbital periods below ~40 minutes. SDSS J1730 shows 
strong helium emission lines with no helium absorption at 
shorter wavelengths; together with the large He I 5875 equiv- 
alent width, this suggests it has an orbital period longer than 
~40 minutes. 

As we push towards fainter targets, removing candi- 
dates that are very unlikely to be AM CVn binaries is im- 
portant in order to ensure that observing the remainder of 
the sample remains practically achievable. Combining new 
constraints provided by GALEX fluxes with our increased 
knowledge of the region of colour space occupied by AM 
CVn binaries, leads to a 43 per cent reduction in the size of 
the sample. This cuts the number of objects still requiring 
observation by 49 per cent, to 275. This should result in a 
corresponding increase in our AM CVn hit-rate, and allow 
the remaining long period AM CVn binaries hidden in the 
SDSS photometric database to be uncovered. This will en- 
able the estimates of the space density to be greatly refined, 
and should reduce the several orders of magnitude uncer- 
tainty in its value. Establishing this will provide constraints 
for common envelope evolution models, allow more accurate 
modelling of the gravitational wave signal and foreground, 
and lead to a better understanding of the stability of mass 
transfer in possible progenitor systems. 
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APPENDIX A: TARGET CATALOGUE 

Our target catalogue of AM CVn candidates selected from 
SDSS DR7, with GALEX UV and Sloan magnitudes, and 
classifications based on our spectroscopic followup, is given 
in Table IA1I The g-ba nd extinction a s sumin g full Galactic 
reddening according to 1 Schlegel et al.l (| 1998T ) is also given. 
It is assumed that at the galactic latitudes of the SDSS, the 
majority of AM CVn binaries will be further f rom the plane 
than t he scaleheight of dust in the Galaxy, see Roelof s~et al.l 
(l2009l ). 
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Table Al. Target catalogue with GALEX UV magnitudes and classifications. Suffixed colons indicate an uncertain classification. Magnitudes given use the AB photometric system. 
This is a sample of the full table, which is available in the online version of this article (see Supporting Information). The full spectroscopic catalogue is available in electronic form at 
the CDS (http://cdsweb.u-strasbg.fr/). 
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